Buenas,

El objetivo de esta pull request es únicamente académico.

He añadido una variable que estaba confusa:
r representa el rango del dominio a mostrar en los gráficos. Estaba confundida con n que a la vez representa el número de aciertos que definen la distribución binomial negativa, BN(n,p), y suponen además un input en el shiny.

Finalmente he modificado los límites del eje vertical del gráfico de la función de probabilidad para que se pueda apreciar la curva de dicha función, aumentando así su carácter didáctico.

Gracias.

saludos

PD. El which.max(aux) lo estaba usando para buscar el valor de la v.a. cuya probabilidad fuera máxima, no obstante, para la finalidad de controlar los límites del eje vertical se puede simplificar a max(aux)